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Abstract 

We study the Hard Core Model on the graphs G obtained from Archimedean tilings i.e. 
configurations in {0,1} G with the nearest neighbor f's forbidden. Our particular aim in 
choosing these graphs is to obtain insight to the geometry of the densest packings in a 
uniform discrete set-up. We establish density bounds, optimal configurations reaching 
them in all cases, and introduce a probabilistic cellular automaton that generates the 
legal configurations. Its rule involves a parameter which can be naturally characterized 
as packing pressure. It can have a critical value but from packing point of view just as 
interesting are the noncritical cases. These phenomena are related to the exponential size 
of the set of densest packings and more specifically whether these packings are maximally 
symmetric, simple laminated or essentially random packings. 
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0. Introduction 



The original Hard Square Model has since its introduction in Statistical Physics surfaced in various 
other fields as well. Apart from modelling the interaction with infinitely strong but zero range 
potential (particles with no interaction except on contact, when they are undeformable), its exclusion 
rule of forbidding two l's in neighboring graph vertices is natural in e.g. communication networks 
and in symbolic dynamics. In these contexts the model is often called independent sets or the golden 
mean subshift. 

In this paper we concentrate on the packing aspects of the model, more specifically on the dense 
packing regime of its two-dimensional version. In short, we are interested in the qualitative properties 
of the densest packings and their formation in a discrete geometric set-up based on uniform tilings. 

First we characterize the geometric and density properties of the optimal packings of l's. After 
this we introduce a probabilistic cellular automaton that is parametrized by packing pressure. There 
may or may not be a dynamic phase transition as this parameter approaches its maximum. Through 
this formulation we will be able to sec the connection between criticality and packing type. In the 
first case there is a phase transition and the densest packings form a finite set of laminated packings. 
The mechanism resulting in the critical behavior is in all our cases the same, a voter rule embedded 
in the hard core rule. The second, non-critical, case corresponds to the set of densest packings 
being exponentially large. This is due to the existence of a local move on the densest packings. In 
this class the optimal configurations are generically random packings. Finally there seems to be a 
qualitatively different "borderline case" corresponding to existence of a non-local move resulting in 
random laminated packings. To further characterize this case as well as some related phenomena 
we also investigate a few cases beyond uniform lattices. 

This work is partly motivated by recent advances in packing in continuous space set-ups (see e.g. 
[CE], [CG-SS]). The overall picture from our 14 lattices compares interestingly with these results. 
Of particular interest is the fact that in R™, n > 1 but "small", the densest packings are believed 
to be lattice or random lattice packings (also in some exceptional dimensions like 24). On the other 
hand for high n as well as in other "big spaces" like the hyperbolic space some type of random 
packings are expected to prevail. In our discrete set-up the random packings surface already in 
planar uniform graphs. 

1.1. Set-up: the Rule and the Archimedean graphs 

Let G be a graph. In all but one of the cases considered in this paper it will be planar. We measure 
distance, d, on G by hop count i.e. by computing the minimal number of edges that need to be 
traversed to move from vertex x to y. The (punctured) neighborhood of nearest neighbors of a vertex 
x is N x = {y G G | d(x, y) = 1}. \N X \ is the vertex degree of x in G. 

Our configurations, which form a subset of {0, 1} G , are defined by the local 

Hard Core Rule: If there is a 1 at the vertex x then on N x the conEguration must be all-0. 



2 



Denote the set of such legal configurations by X^f. 

Most earlier studies on the hard core model have concentrated on the square lattice or more generally 
on Z d . Additional results exist on the triangular lattice, on trees and a few other set ups, see e.g. [B], 
[GK], [BW]. While the model on a graph is the most general one, it is not the natural set-up when 
one is interested in the packing aspect. Graph as a purely topological object needs to be augmented 
with a metric structure so that the nearest neighbors are in some sense also geometric neighbors. 

Z 2 , T (triangular lattice) and H (honeycomb lattice) are the vertex and edge sets of the three 
regular tilings of the plane. Their most immediate generalizations arise from the uniform or 
Archimedean tilings. While the regular tilings are constructed by tiling the plane with a single 
type of regular convex polygon, for uniform tilings a mixture of the different types of these polygons 
is allowed. Turns out that there arc 11 such tilings of the plane. They share the common property 
that at each vertex the tiling looks up to rotation identical and the nearest neighbors (one edge 
away) are all at at the same Euclidean distance away which we take to be the unit. The code names 
of the tilings are obtained by circling a vertex and recording the n-gons along the way. Hence Z 2 
is 4 4 , the triangular lattice is 3 6 etc. For more on properties of the these tilings see [GS]. From 
now on we shall view the vertex and edge sets of the Archimedean tilings as geometric graphs, the 
Archimedean graphs/lattices (most of them are not mathematical lattices but this terminology 
is in use nevertheless). Figure 1. illustrates their set. 

1.2. Densest packings 
1.2.1. Maximum density 

Let us first investigate the density and geometric properties of the densest packings on the uniform 
lattices when the Hard Core Rule is imposed. 

The density of l's in a configuration c is 



where the limit must exits and agree for all sequences {D n } retaining their two dimensional shape 
as n — > oo (D n is a domain of size n lattice points on G). This number, when existing, is always 
between and 1. 

For some of our graphs it is immediate which sets support the densest packings (e.g. for Z 2 the 
checkerboard), for others it is less obvious. We first argue the optimal densities and present samples 
of corresponding configurations. 
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(2) 

Consider a vertex x € G and the set of its second neighbors i.e. vertices two hops away, 

c(x) = 1 then by the rule on N x the configuration is all-0. On the set N x there may be both O's 

and l's. For the purposes of the upper bound argument we will ignore the configuration outside 

(2) (2\ 

{x} U N x U Nx ■ For y £ N x let n y be the number of y's nearest neighbor l's inside Nx and n be 
its average: 

(1.1) n = — \- - max V" n v . 

(2) 

Note that n|-/V x | always majors max|{z e iV^ | c(z) = 1}| which is the discrete analog of the 
kissing number of the sphere packing (see e.g. [CE]). 

Proposition 1.2.1.1: On an Archimedean graph of degree d and second neighbor mean occupancy 
n the global conhguration density of 1 's is bounded from above by 

(1-2) P=-r A ir ■ 



' n+l 

(2) 

Proof: For the 1 at x there are d O's neighbored by at most dn l's in Nx ■ Hence the average 
number of l's each is shared by is at most n+l. Equivalently, on the average for each 1 there 
are at least d/(n + 1) O's. On an Archimedean graph d is constant and since all neighborhoods are 
rotations of each other n is constant, too. Hence the bound holds at every vertex on the graph and 
the upper bound follows. I 

Remarks: 1. The result is actually more general than our set-up and also applies e.g. to homoge- 
neous trees: n = d — 1, hence p = 1/2 always. 

2. Despite its simplicity the bound is often tight. For example for the square lattice we have d = 4, 
n y = 3 hence n = 3 and p = 1/2. Indeed it is tight for the first five and last two Archimedean graphs 
as listed in Table I. 

For the remaining Archimedean cases the density bound is argued as follows. 

Proposition 1.2.1.2: For (3 2 .4.3.4), (3 4 .6), (3 3 .4 2 ) and (3.4.6.4) the maximal density is 1/3. 

Proof: The infinite limit density bounds follow from the following arguments for configurations on 
a finite torus. 

(3.4.6.4) : the graph is made of disjoint unit triangles (all unit edges) and since each can support at 
most one 1, the density is at most 1/3. 

(3 2 .4.3.4): each vertex is shared by three unit triangles, hence in the configuration the number of 
vertices is the same as the number of triangles. Since each triangle can carry at most one 1 the given 
bound follows. 

(3 4 .6): each vertex is shared by four unit triangles, hence the number of triangles is 4/3 of the 
number of sites. Therefore the number of l's is bound by 1/4 of the number of triangles which 
implies the bound 1/3. 



(3 3 .4 2 ): in each horizontal strip of triangles every 1 forces two O's to its right implying the bound. I 
Examples of packings reaching the maximal density are shown in Figure 1 (l's are solid dots). 

1.2.2. Qualitative features 

Table I. lists some of the qualitative and quantitative packing properties of the Archimedean graphs 
(the first 11 lines) and three related set-ups that we analyze later. For Z 2 and T existing results 
from the literature have been utilized but essentially all other data presented in Table I is established 
in this paper. 

The first three lines correspond to the regular lattices. For example the first line indicates 
that the optimal packing on the square lattice (4 4 ) is supported by one of two square lattices (with 
lattice spacing tilted by 7r/4, hence diamond tile) which as sublattices each have density 1/2. As 
a packing it is simply what we would get by stacking the densest packings on one lower dimension, 
Z, therefore it is denoted by L for a laminated packing. 

The first seven cases are similar, call them rigid: for each graph there is a (small) finite number 
of subgraphs all isometric to each other which support the densest packing and these arrangements 
of l's cannot be deformed in any fashion. Topological equivalence means that there is a home- 
omorphism between the subgraph and the stated Archimedean graph (in the terminology of [GS] 
the two graphs are of the same topological type) . In all our cases the homeomorphism is actually 
close to the identity map so the actual subgraph "looks like" its preimage. See also Figure 1. for 
the geometry of the subgraphs involved (l's are indicated by the solid dots). 

For (3 3 .4 2 ) the densest packing is an infinite family of random laminated/slide packings. 
If every third entry in the 2-rows made of the triangles (framed in Fig. 1) carries a 1 they can 
be stacked with off-sets ±1/2. Hence in a square domain ofnxn vertices there are approximately 
e h<1) ™, — |log2 different packings all with p = 1/3 (the superindex of h refers to this being a 
1-dimcnsional entropy). 

The remaining three cases (3.4.6.4), (3.12 2 ) and (3.6.3.6) (Kagome lattice) are yet qualitatively 
different from those above. In all three the densest packings are essentially random packings. 
They all have p = 1/3 packings that allow local moves: local exchanges of O's and l's between 
neighboring lattice sites (which preserve the legality of the configuration). Because of the locality 
of this action the number of packings in a domain ofnxn vertices grows like e' i<2> " 2 , > 
depending on the graph. For (3.4.6.4) and (3.12 2 ) see Fig. 1 for the details of these exchanges (the 
arrows) . 

Subsequently we call maximally symmetric densest configurations optimal packings. For lami- 
nated packings these are all there are, for the other two types they are rare specialties. 

Having identified a densest packing subgraph for an Archimedean graph and viewing it as a subset 
of R 2 wc can define its Voronoi tessellation. Every 1 is the center of a convex 1-tile which is unique 
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up to rotation. If the subgraph is in fact one of the Archimedean graphs this tile has the O's of 
the set N x along its edges. An identical tile with the symbol 1 at the center replaced by is a 
O-tile. In Table I we have identified the types of these tiles and some of them are illustrated in 
Figure 1 (in bold line) . Archimedean graphs for which this construction yields the simplest tiles are 
(4 4 ), (6 3 ), (3 6 ) and (3.6.3.6). Note that considering the packing of these 1-tiles on R 2 instead of l's 
on G does not introduce an extra condition. Because of the Hard Core Rule there cannot be 1-tilc 
overlaps. If a tile arrangement in plane is in fact a tiling i.e. a one-fold cover of the plane then it is 
necessarily the densest packing of l's. By the arguments above the rigidity properties of the densest 
packings translate into the rigidity of these tilings. 



2.1. Dynamical packing with a PC A 

We now introduce a dynamics that enables us to generate increasingly dense packings in all our 
set-ups. This yields further insight into their nature and in particular will clarify the nature of the 
three qualitatively different packing classes identified in the previous section. 

Compatible with the Hard Core Rule the configurations can be changed locally according to 
the following local 

Update: If there is at least one 1 on N x then the symbol update at x is 0, otherwise it is 1 with 
probability p independently of the updates outside N x U {x}. Denote this random local map from 
N x to {0, 1} by f p . 

Instead of the probability p a related quantity called the activity (fugacity, z — p/(l — p)) can be 
used. Sometimes it is taken to be sublattice dependent corresponding to sublattice preference (akin 
to an external field in the Ising model). We do not consider this variation; our rule is uniform over 
the lattice. Its a simple but important result that the update, when nontrivial, is an irreducible 
action on the configurations: 

Proposition: For < p < 1 the Update can transform any legal hnite configuration to any other 
such configuration with a finite number of steps. 

Proof: Given any legal finite configuration, there is a positive probability that it will be transformed 
by the update in one step to the all-O-configuration. Since this operation is reversible for the given p 
any two such configurations can thus be transformed to each other with positive probability in two 
steps. I 

An efficient way of generating all legal configurations from any given configuration can be given once 
the subgraphs carrying the optimal configurations have been identified (as in Table I): 

1. Take any initial configuration on G (legal or not). 

2. Pick any optimal subgraph and apply the Hard Core Update at neighborhoods centered at its 
vertices to generate an updated configuration on the subgraph. 
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3. Cycle through all the optimal subgraphs in turn using the rule in 2. Once this is done all vertices 

of the graph have been updated. 
4- Combine the subgraph configurations into a global configuration on G (which will now be legal). 

The global map F p : X^f — ► X'^f. defined by steps 2-4 is our Probabilistic Cellular Automaton 
(PC A). Iterating it will relax the initial configuration to a Gibbs state on X^f (a measure of maximal 
entropy on configurations compatible with the Hard Core Rule). 

Example: On T the three sublattices are shifts and tilts of V3T, call them dot, circle and ring 
sublattices (optimal sublattices, identified this way in Figure 1). If we take e.g. all-1 initial state on 
T (illegal) and choose first to update the dot sublattice, it will become all-0 independent of p (since 
circle and ring are still all-1 and N x consists of them for each x in dot sublattice). In step 3 dot and 
circle will in turn update the ring sublattice the same way to all-0 (due to l's on circle). This is 
followed by circle being updated by dot and ring according to the probability p. After the first full 
cycle through the three sublattices the subconfigurations combine to a legal configuration. 

Among our graphs either two or three sublattices are involved in the iteration except for Z 2 M, which 
cycles through four sublattices. 

By the Proposition the action of the non-trivial PCA is irreducible on the set of Hard Core con- 
figurations. Hence in particular this algorithm can be used to generate all configurations from one 
sample (in a finite number of steps for a finite lattice). 

It is natural to view the update probability p as the packing pressure: the PCA with higher 
p is just simply squeezing more l's into the available slots. In terms of the auxiliary tilings just 
introduced we fill in new 1-tiles centered on the subgraph being updated if the rest of the 1- tiles 
allow i.e. if no overlaps result. 

2.2. Criticality and noncriticality in packing 

The PCA above has been coded for our set of lattices and we will present here the results from 
the relaxation study. Apart from the square and triangular lattice cases which can be found in the 
literature the subsequent data is from our runs. 

The aim here has not been great numerical accuracy but rather in verifying the critical behavior 
with a reasonable estimate for the critical parameter p. We have used the following sequential 
procedure. First we relaxed on a small (about 50 x 50 sites) toral lattice to see if p c < 1. If so 
we then obtained a better estimate for p c on a larger toral lattice (at least 200 x 200). Finally we 
employed "worst case" initial conditions. In the subcritical regime we relaxed initial conditions filling 
completely one of the sublattices to see if the sublattice densities converged to a common value. In 
the supercritical case we used Bernoulli initial distribution with density of l's approximately p(p c ) 
and observed if the densities on the sublattices asymptotically diverged. This procedure was done 
with an increasing/decreasing sequence of p's in the sub/supercritical case respectively to obtain a 
bracket for p c . 
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2.2.1. Regular and other lattices with laminated packings 



The model on Z 2 is sometimes referred to by the name the Hard Square Gas. The squares are 
actually diamonds of side length y/2 that we have identified above, in Table I and in Figure 2. The 
critical value for p, p c w 0.79 was computed to high accuracy in [BET] and the Dobrushin-Shlosman 
criterion for phase uniqueness was applied in [RS] to provide a rigorous lower bound. Below the 
critical probability the model is ergodic, above the dynamics preserves two Gibbs measures assigning 
different probabilities on the two sublattices. 




On the honeycomb lattice H the dynamical behavior is similar to the square lattice case (for critical- 
ity proof see [R]). Here the geometry is as indicated in Figure 1: hexagonal 1-tiles on two triangular 
lattices. In the supercritical case the densities of 1-tiles are unequal on the two sublattices whereas in 
the subcritical case they converge to a p(p). Figure 3. shows a typical relaxation of these densities in 
the subcritical regime. The mirror-effect in the graphs is due to the fact that a high packing density 
(p w 0.384) has been attained and preserved in spite of the fluctuations around the equilibrium level. 
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Figure 3. Relaxation on the honeycomb lattice (p = 0.859, 200 x 200 torus, 5.000 iterates.) 

On the triangular lattice T the model corresponds to the Hard Hexagon Model of [B]. It is known 
that beyond p c s=s 0.92 there is a three way sublattice split into identical but thinner triangular 
lattices supporting the hexagons. Our runs indicate that the highest density of a packing that does 
not yet reveal the sublattice parity is approximately 0.26. 

Our simulations on 200 x 200 site (sometimes bigger) toral lattices established that on the 
lattices (4.8 2 ) and (4.6.12) the two sublattice alternatives carrying the densest packings emerge 
above critical values p c 0.9 and 0.92 respectively. In the case of lattices (3 2 .4.3.4) and (3 4 .6) three 
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way splits take place approximately at the critical probabilities 0.99 and 0.97. The corresponding 
critical densities are listed in Table I. 



2.2.2. Embedded voter rule 



The critical behavior in the cases above seems to result from a common underlying structure. In 
short in these models one can identify a sublattice on which the square of the update is behaving 
essentially like a slightly asymmetric voter rule. 

Suppose that we have configurations from {0, 1} G and a finite neighborhood N. The classical 
majority voter rule has the following basic properties (see e.g. [L]). 

• If N has n s copies of s then the update to this symbol takes place with a non-decreasing 
probability p(n s ). If n s > \N\/2 then p(n s ) > n s /|iV| i.e. the update favors the majority. The 
update is symmetric in the symbols and if \N\ is even then in case of a tie in N the symbol is 
updated without a bias i.e. w.p. 1/2. 

• The updates in disjoint neighborhoods are independent. 

Intuitively these mean that the interpolated graph of the function p(n) is "/ -shaped" , symmetric 
with respect to (|7V|/2, 1/2) and hits the corners (0,0) and (\N\, 1). 

Consider the map / 2 on the honeycomb lattice. It is a random map from a 7-tuple to the set {0, 1}. 
The 7-tuple is the dotted subset on the v^T sublattice as illustrated in Figure 1. Similarly / 2 e.g. 
for Z 2 is supported by a 3 x 3 diamond on a \/2Z 2 sublattice. 

Now take p = 1 and plot the fraction of neighborhoods that yield update 1 under j\ as a 
function of the number of l's in the neighborhood (the 7-tuple or 9- tuple identified above). These 
have been plotted in Figure 4a and 4b, the leftmost graphs in each. The graphs clearly have the 
same qualitative character as the voter curves except that they favor the update 1 even stronger 
resulting in slight asymmetry. 

In the triangular lattice case the set-up is slightly more delicate. In Figure 1 the dots, circles 
and rings in a hexagonal array identify the 19-tuple from which the update to the center dot is 
determined. In this case we consider the restriction / 2 of / 2 to the dot sublattice. This map from a 
7-tuple is again a voter- type map as seen in Figure 4a, right. The other sublattices work identically. 
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Figure 4a, b. Voter curves for H, T, Z 2 and Z 2 M (2) graphs (left to right) 
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The key idea here is that the voter rule is known to be critical ([L]) i.e. depending on the steepness 
of the graph of / 2 around | AT|/2 there is p c < 1 s.t. above it the dynamics becomes nonergodic. 

It is easy to see that for / 2 , p w 1, the update probability closely approximates the three curves 
above. Each one of them is for high p a rule that favors the l's to take over one of the sublattices 
thereby blocking the l's from the other sublattices. Hence the criticality in these models seems to 
coincide with the criticality of the majority voter rule. 

2.2.3. Slide packings 

Among the Archimedean graphs our model on (3 3 .4 2 ) is exceptional. Table I. indicates that its 
high density characteristic property is the randomness in the lamination i.e. stacking of 3Z densest 
packings. We now investigate this phenomenon by going into a slightly richer lattice set-up. 

Let us reconsider the Z 2 lattice with the Moore neighborhood i.e. with neighborhood set extended 
to the nearest eight neighbors in the Euclidean distance on R 2 . The neighborhood graph is not 
planar anymore since the graph edges of the diagonals do not intersect. 

In the hard core model on this graph the 1-tiles are 2x2 squares. There are four maximally 
symmetric tilings each supported by a copy of 2 Z 2 . They are optimal packings but not the only 
ones. In such a packing we can pick any column/row with a 1 in it and shift the entire column/row 
by ±1 lattice sites. This is exactly the kind of behavior that distinguishes the model on (3 3 .4 2 ), too. 

But Z 2 M is a bit more subtle. After shifting one column we can still slide in similar fashion any 
other column but none of the rows. This obviously generates an infinite number of configurations 
that are all densest possible. The one-dimensional entropy h\ = linin^oo iln2 n / 2 = 5 In 2 is the 
exponential growth rate of the number of these configurations from infinite strips of width/height n 
lattice sites. 

In terms of packing the optimal ones are still all laminated packings but somewhat nontrivially 
so. A pair of sublattices that are exactly one unit lattice shift (vertical or horizontal) from each 
other can together support a perfect tiling whereas a pair of sublattices that are one horizontal and 
one vertical unit shift from each other can not. Because a random mixing of the two sublattices 
is allowed (as in the counting above) now essentially all of the optimal packings are disordered 
laminated packings. Some of their (non)ergodic properties were analyzed in [E]. 

From the above we can expect that the equilibrium of the PCA evolutions when p "f 1 is more 
complex than in the previous cases. This is indeed the case and in particular understanding the 
critical behavior becomes more challenging. There are an infinite number of ground states but in 
fact the only critical transition is between the two possible "orientations" of the ground state. This 
has caused significant difficulty in earlier studies (like [BN] , where restriction to an infinite cylinder 
broke the symmetry thereby obscuring the criticality analysis. For a more recent analysis of this 
model see e.g. [LC]). 
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In our simulations there seemed to be clear phase segregation above the approximate threshold 
p c = 0.98. There a pair of sublattices take over i.e. the tiles increasingly concentrate on them 
deserting the two others. The pairing is forced by the structure of the densest packings above. 1-tile 
configurations on such a sublattice pair combined can together achieve arbitrarily good cover of R 2 
and indeed tile it. 

Figure 5b illustrates the near equilibrium configurations on the four sublattices in the super- 
critical regime (p w 0.992). The four possible sublattice matchings are the pairs next to each other 
in the same row or column in the 2x2 array shown. Clearly here the two dominant ones are the 
ones on top row which superimposed to each other approximate fairly well the p = 1 limit of a row 
shifted disordered tiling. 



One should expect the voter mechanism to be more complex in the context of the Z 2 M lattice. As 
we have seen, there any one of the four sublattices can support only a single densest packing (one of 
the maximally symmetric ones), whereas there are four pairs of sublattices that can each support an 
infinite number of densest packings. The criticality emerges from the system's decision on whether 
the horizontal or vertical slides are allowed in these packings. 
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Figure 5a, b. Z 2 lattice with Moore neighborhood squared and a supercritical state. 
(p = 0.992, 200 x 200 torus, 10.000 iterates from S(l/8).) 



Figure 5a indicates the 25-tuple from which the center site is updated under / 2 . Form now doublets 
from (horizontal) nearest neighbor dots (one black and one white). / 2 induces a unique update map 
on these, denote it by / 2 Q . The entries of the doublet (»o) as well as the update can be of the form 
(0, 0), (0, 1) or (1,0). / 2 is again a map from 9-tuplcs (of doubles) to {0, l} 2 and moreover by the 
Hard Core rule among the 6x3 lattice points on which / 2 Q depends there are at most nine l's. 

The middle curve in Figure 4b gives the fraction of 3 x 3 non-overlapping doublet neighborhoods 
with k l's in them that update to a doublet with 1 in it. This map is identical to the three maps 
considered above, only now we map horizontal doublets instead of singletons (individual lattice site 



12 



values). The reason that one must consider / 2 Q instead of / 2 is that this map is independent of the 
horizontal slides in the configuration. It's graph indicates a weak majority voter character. 

The rightmost curve in Figure 4b is what one gets by restricting f\ to one sublattice alone. The 
map clearly fails to be majority voter-type. If one defines doublets by combining sublattices that 
are one horizontal and one vertical shift from each other the resulting map is not voter-like either. 
This is expected: this arrangement of 1-tiles leaves gaps in the packing and cannot prevail under 
high packing pressure. 

The above seems to support observed critical phenomena in the models. It shows that a majority 
voter mechanism operates at the level of the optimal sublattices. There are caveats though. As 
seen from the graphs the f\ maps are not quite symmetric as they should be in the ideal voter 
rule. However this slight asymmetry actually favors the phase segregation. Spurious combinatorial 
phenomena may show up. There may e.g. be "instability " of the following kind: Consider f\ on 
Z 2 and suppose the center entry in the 3x3 diamond is 1 and the rest of the entries are O's. This 
neighborhood yields the update 1. But if the entry 1 is anywhere else in the neighborhood while 
the others are still O's the update will be 0. Whether details of this kind can influence the global 
relaxation of the PCA for p < 1 will require further study. 

2.2.4. Noncritical packing 

The third and perhaps the most interesting class of packings on a lattice are the 2-d random packings. 
This class seems to coincide with the noncriticality of the PCA, a connection of definite interest. 
The analysis here concentrates on the Kagome case but we believe that the same principles apply 
verbatim to the three cases among the Archimedean set and beyond. 

Kagome lattice (K) is a 1/2-thinning of the triangular lattice obtained by removing every other 
lattice line in all of the three lattice directions. It has appeared in the context of a number of 
statistical physics models ([B]) and it turns out that here, too, it yields some interesting insight in 
to the nature of the high density packings on lattices. 

Each lattice point has four nearest neighbors in a "hour glass" formation. If we would choose 
the 1-tile to be a hexagon with N x among its vertices these hexagons distributed according to the 
l's from a legal configuration could never overlap, but they would not form a tiling either. However 
if one extends the hexagon to a rhombus with long axis aligned to the vertical axis of the hour glass, 
then these 1-tiles cannot ever overlap and moreover will form a tiling for the densest packing (a 
rhombus is two equilateral triangles glued together at their base). This construction immediately 
implies that the maximum density of l's on the lattice is 1/3. 

The lattice divides in a natural way up to three identical sublattices, each of which is again a 
Kagome lattice scaled by factor and rotated appropriately as indicated in Figure 1 (the dots, 
circles and rings identify the three sublattices). Each of these sublattices supports a rhombus tiling. 
The arrangement on the left, three rhombi forming a hexagon, can be viewed as a double covering 
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set since with it one can cover the plane with one rhombus overlap between neighboring hexagons. 
All the rhombi in this cover are supported by the same sublattice (the dots). These three covers 
correspond to the most symmetric of the densest packings on the Kagome lattice. 

The three-tuple of rhombi forming a hexagon can be rotated by tt. This flip is an example of a 
local move. In the new formation the rhombi are centered on the ring sublattice. Hence in a domain 
containing N non-overlapping hexagons there are at least 2 N different rhombus tilings. Indeed the 
exponential growth rate of rhombus tilings (entropy per tile) can be shown to be 



(which cannot be evaluated in closed form but is approximately 0.32306, see [DMB], [W]). The set 
of rhombus tilings is closed under flips in the sense that any two tilings can be reached from each 
other through a finite sequence of flips (for a most general result in this direction see [P]). 

The case at hand is very different from those in the previous sections including the Z 2 M. In the 



whereas on Kagome this is not the case: if weighted uniformly almost all of them are disordered 
without any periodic structures left in them. In other words in this class the densest packings are 
generically random packings. 

The relaxation dynamics of the Hard Core PCA defined on Kagome lattice is unlike in the cases 
above. There is no majority voter mechanism to favor any sublattice and there is no critical value 
for the packing pressure p. In our simulation runs (about 200 x 200 site tori) the relaxation appeared 
to be exponentially fast for all p and the limiting configuration was evenly distributed on the three 
optimal sublattices. This is to be expected: the equilibrium measure should be the unique measure 
of maximal entropy since this is unique even at the p = 1 limit (where it is the uniform measure on 
the optimal packings, additionally also invariant under the flips). 

Hard core on (3.4.6.4) and (3.12 2 ) are qualitatively similar to the above in high density packing. 
Figure 1 indicates densest packings with p = 1/3. Both admit 2d random perturbations with the 
indicated local moves (arrows). We have computed the entropy bounds in Table I for (3.4.6.4) from 
the independent choices on the period parallelogram (three inside the shaded indicated domain) and 
for (3.12 2 ) from the rotations in density 1/2 12-gons in the period parallelogram. 

2.2.5. Extensions: Union Jack and Quilt lattices 

Union Jack (UJ) is derived from the unit square lattice by adding the ascending and descending 
diagonals. It can also be obtained as a dual of the Archimedean lattice (4.8 2 ) (i.e. it corresponds 
to the Laves tiling [4.8 2 ], see [GS]). Quilt, Q, for which we do not know a more standard name, is 
the thinning of the Union Jack when one removes every other ascending and descending diagonal. 
These lattices are illustrated in Fig. 6a and b. 




last one the number of densest packings in a square of n 2 vertices is asymptotically e' l<1) " whereas in 
the Kagome case it is e' l<2) ™ 2 . The reason for this is that on Z 2 M all the packings are still laminated 
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Figure 6a, b, c. Union Jack and Quilt graphs with tiles. Approximate critical curves. 



The interest in these lattices stems from the fact that they refine Z 2 in two new ways and in the 
process bring along a novel type of 1-tile. In the case of the Union Jack lattice it is a unit square 
(hence having half the area of the diamonds) and in Quilt another diamond. 

Our update rules will now have two parameters, the probabilities P4 and either p s or p 6 . The 
index refers to the degree of x, squares having four and diamonds eight/six neighboring sites on 
their boundaries in this lattice in Union Jack/Quilt respectively. 

The model with p 4 = is the Hard Square Model. The critical behavior on both lattice set- 
ups extends to the pn > regime. In simulations the phase diagram appears approximately like 
in Figure 6c, the curves separating the sub- and supercritical phases. The curves should be above 
the p c w 0.79 horizontal line, increasing, reaching the northeast corner and be ordered as shown 
for the two models. Our simulations (on 100 x 100 torus) support these claims, which can also be 
heuristically argued but we do not have a rigorous argument. 

From the packing point of view Union Jack lattice is somewhat ambiguous since the squares and 
diamonds do not have same area. In the high density limit p±, p$ | 1 all-l-squares arrangement should 
prevail. Hence among the optimal packings there should be a unique densest one (off pg = 1). This 
behavior disappears in the Quilt lattice - the squares are there replaced by diamonds of degree four 
which are of equal area to the degree six diamonds. The optimal packings for the Quilt lattice are 
clearly the three sublattice tilings. No other perfect tilings exist in this set up. 

The models on the Union Jack and Quilt lattices are for small p 4 just variants of the model on Z 2 . 
Consider again the 3x3 diamond neighborhood on Z 2 . If in each of its 2 x 2 subneighborhoods 
there is at least one symbol 1 then on can compute that whereas / 2 yields update 1 on Z 2 w.p. p, 
on Union Jack the update with / 2 p4 p ^ will be 1 at least w.p. p&(l — p 4 ) 4 and on Quilt at least w.p. 
Pq(1 — pi) 2 . Hence the update probability on both of these lattices is near to the square lattice case. 

3. Conclusions 

In this paper we have studied the densest packings of l's on certain planar lattices with the ex- 
clusion implemented by the nearest neighbor Hard Core Rule. The choice of uniform lattices was 
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paramount: we wanted to investigate the most natural discrete analog of the densest packing of 
disks/spheres/hyperballs in the Euclidean world. 

As we have seen a rather striking picture emerges; either there are only a few densest packings 
or there is an exponential number of them (in the size of the domain) . This split is a consequence 
of the densest packing being either rigid or loose (i.e. the ground state having zero or positive Id 
or 2d residual entropy as for example in the Ice model of Statistical Mechanics) . This sparseness or 
richness of the set of the densest configurations in turn is a direct consequence of the existence of a 
local move. The lack of a local move implies long range order whereas the existence of a local move 
enables one to join densest finite patches into a densest global configuration. 

The packing dynamics can be critical and the packing types are intimately related to this criti- 
cality/noncriticality in the generating PCA. In all our cases the lattices with rigid ground states (no 
local move) are approached via a critical transition when the packing pressure increases. Conversely 
when the set of ground states is rich enough, they can be piecewise glued together to a global one 
and a critical transition is not imposed on the model in the quenching of the PCA. 

The processes involved here are quite subtle which is perhaps best indicated by the existence 
of the borderline critical case of slide packings. These are likely not to be combinatorial oddities - 
indeed they resemble a great deal the well known Barlow packing, the optimal packing of 3-d hard 
spheres (see [CG-SS]). In this context there is no local move but there is an infinite slide move. Hence 
there is limited 1-d long range order and apparently a phase transition in the packing pressure that 
goes with it. There has been attempts to clarify this using e.g. the Pirogov-Sinai theory to no avail 
(private communication) . To fully understand the nature of the transitions in these models more 
studies are needed. 
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Table I. Packing properties of the Archimedean and other graphs. The first three lines cor- 
respond to the regular graphs, the first 11 to the Archimedean ones and last three are non- 
Archimedean. 

• Optimal subgraphs: the named ones carry the most symmetric densest packings. ~ means 
topological equivalence. 

• Multiplicity: number of densest packings. Finite number means that all optimal ones are 
isometric to the most symmetric one, infinite that there are random perturbations. 

• Tile: n-gon is not regular, whereas diamond, square and hexagon are. Rhombus forms from 
two equilateral triangles glued together at base. For further details see Section 2.1. 

• Type: L = laminated packing, RL = random laminated packing i.e. Id slide packing, R = 2d 
random packing. 

• p: density of the optimal packings in the full graph. 

• h^: the one- and two-dimensional residual entropies of the densest packings. 



